New Observables in Break-up Reactions of Three-body Halos. 



E. Garrido, D.V. Fedorov and A.S. Jensen 
Institute of Physics and Astronomy, 
Aarhus University, DK-8000 Aarhus C, Denmark 
(February 5, 2008) 



o 

=5 



X 



C3 



Momentum distributions of particles from fast three-body halo fragmentation reactions with po- 
larized beams and outgoing fragments are analyzed by use of the sudden approximation. The final 
state interaction between the two non-disturbed particles is considered. We first give a general 
| and detailed description of the method. We introduce observable quantities that emphasize the 

. two-body correlations in the initial three-body structure as well as in the final two-body system. 

Q\ ' Using neutron removal from 11 Li as an example we investigate the dependence of the polarization 

t-H t observables on the properties of the low-lying resonances and virtual states of the 10 Li subsystem, 

^jjy These observables are very sensitive to the I > 0-waves in the wave functions, and they provide 

detailed information on the resonance structure of the neutron-core subsystem. 

<\ 

I PACS: 25.60.Gc, 25.60.-t, 24.70. +s, 21.45. +v 



> ' I. INTRODUCTION 

o : 

Halo nuclei have been extensively investigated in the latest years both from the theoretical and experimental point 
qq ' of view. For a review see for instance and references therein. 

The general properties of light nuclear halos have succesfully been revealed by means of few-body models, dividing 
VO ■ the de grees of freedom into the approximately frozen (core) and active (halo). Special attention has been paid to 
three-body Borromean bound systems [0-p| , where none of the two-body subsystems is bound. The most prominant 
examples are 6 He ( 4 He+n+n) and 11 Li (^Li+n+n) both extensively discussed in a general theoretical framework [0. 

The most detailed source of experimental information about the structure of these nuclei is the momentum distribu- 
tions of the particles resulting from the fragmentation reactions |t]-[u| . For one particle removal reaction with a high 
energy beam we can consider that the removed particle in the halo projectile is instantaneously captured by the target 
without disturbing the remaining two. This reaction picture is known as sudden approximation, and has been proved 
to be fairly accurate in the description of such halo fragmentation reactions Jl4|-p9[. For a correct interpretation 
of the experimental data it is essential to include the interaction between the two non-disturbed particles after the 
fragmentation, especially when low- lying resonances are present [ p~T] , p~5|Jl~^ ] . Detailed and consistent incorporation of 
this Final State Interaction (FSI) into the model is necessary as seen by its influence on the momentum distributions 
p9| . Neutron momentum distributions are very much affected by the final state interaction, but also core momentum 
distributions are significantly influenced. 

A polarized beam interacting with a target is not equally sensitive to all the components of the wave function. It is 
then possible by means of this tool to observe signals of the less dominant components, and in this way to investigate 
finer details of the nuclear structure of the projectile. The purpose of this paper is to investigate the effects of polarized 
beams in fragmentation reactions of fast three-body borromean halo nuclei. The case where the spin projection of 
one of the particles is measured in the final state will also be considered. We introduce several useful observables and 
show the kind of information we can extract from them. In particular we consider neutron removal from the polarized 
11 Li projectile. The resulting neutron momentum distributions are sensitive to the amount of p-state admixture in 
the neutron-core relative state. 

The correct description of the structure and reactions of three-body halo nuclei is an indispensable prerequisite 
for the present investigation. These general procedures are discussed in previous papers P,p0|-p2| which we shall use 
extensively. This paper is organized as follows: In section II we formulate the theoretical model and the method 
used to describe the fragmentation reaction including definition of the polarization observables. The numerical results 
obtained for 11 Li are reported in section III and in section IV we give a summary and the conclusions. We collect the 
rather complicated analytic expressions of the differential momentum distributions for different reactions in appendices 
A and B. 



1 



II. METHOD 



The fragmentation reaction is described in the sudden approximation, which assumes that one of the three particles 
in the projectile is instantaneously removed by the target without disturbing the remaining two. This approximation 
is valid for a high-energy three-body halo beam interacting with a small target through short range interactions. 
Therefore we assume a light target where the Coulomb dissociation process only contributes marginally. 

The transition matrix in the sudden approximation is given as 

M = ($|*) , (1) 

where $ is the final state wave function after break-up and ^ is the initial wave function of the three-body projectile. 
We now must specify these wave functions, compute the matrix element and integrate \M\ 2 over all the non observed 
variables. 



A. Wave functions 



The initial wave function describes the bound three-body halo system whereas the final state wave function describes 
the two-body continuum state specified by the experimental conditions. 



1. Initial wave function 

The projectile three-body wave function ^ with inclusion of the proper quantum numbers and variables is denoted 
by i$/ JM (x, y), where J and M are the total spin and its projection. The coordinates x and y are the usual Jacobi 
coordinates pQ ], where x is drawn between the two particles surviving after the fragmentation. The wave function 
is expanded in terms of a complete set of hyperspherical harmonics Ygjt (a,n, x ,£l y ), where the quantum number K 
usually is called the hypermoment, l x and l y are the orbital angular momenta associated with x and y, and L is 
resulting from the coupling of these angular momenta. The variables p, a, Q x , and fl y are the hyperspherical variables 
obtained from the Jacobi coordinates x and y B. Then ty JM (x,y) is written as 

* JM (x,y) = 

E CnKi x t y LsM [y&(«.fix,n»)®xf.. y Y M ' (2) 



p 



where the spin function Xs mSy i s obtained by coupling the spins of two of the particles to s x which in turn is coupled to 
the spin s y of the third particle resulting in the total spin S. The radial functions f n (p) and the expansion coefficients 
CnKtxiyLsxSip) are calculated numerically by solving the Faddeev equations in coordinate space p0[ . Details about 
this procedure can be found in M. 



2. Final state wave function 



The final state wave function <!> in the plane wave approximation in the center of mass frame of the projectile is 
given by 

$ = ^k.R^iq.r^^ 5 (3) 

where xV f° r i=l,2 and Xsy are t ne spin wave functions respectively for the two undisturbed particles and the third 
removed particle. Here s and a are the total spin and its projection on an axis. Furthermore k and q are the total 
and relative momentum of the two remaining particles in the final state, r is the distance between the two remaining 
particles, and R is the distance between the center of mass of this two-body system and the removed particle. 

The description in terms of plane waves in Eq. (^) implies that all the interactions between the particles are neglected 
in the final state wave function. However, in our approximation only one of the particles is removed from the other two, 
and the interaction between the non-disturbed particles can not be neglected. In a recent work ]l9t we investigated 
the influence of this final state interaction (FSI) on the momentum distributions, and observed that they are crucial 
for reproduction of the experimental momentum distributions. 



2 



In Eq.(||) introduction of the final state interaction amounts to substitution of the plane wave, e l( *' r X°lX%% i de- 
scribing the two-body system in the final state by the appropriate distorted two-body wave function w. When the 
two-body interaction do not mix the two-body states with different spin s x and relative orbital angular momentum 
£ x , we can expand this two-body wave function (w) in partial waves as in p3| . These assumptions are usually strictly 
valid. The only exception arises from the tensor interaction. However the resulting mixing is often very small and 
therefore insignificant in the present context. We then have 

w aia2 [k x ,x) = ^2(s 1 a 1 a 2 a 2 \a x a a! )w'»' r "(k x ,x), (4) 

with 

w ^( kx , x ) = ] H±- y utJk x ,x)y< xSx (n x ) 

x y {^xmt x s x a x \j x m x )i ix Y tio1n ^{VL kx ) , (5) 

where fc^ is the momentum associated with x, a x is the spin projection of s x , u is the radial and y the angular 
distorted wave function. The angles Q x and fifc x define the direction of x and k x , respectively. The dependence on 
the spins s± and s 2 of the two particles surviving after the fragmentation reaction has been omitted for notational 
simplicity. 

The radial wave functions u J f s {k x ,x) are obtained numerically by solving the Schrodinger equation with the 
appropriate two-body potential V(x) 

^tsS^x) + - -jpV&Jx) - xKx x2 > ) u J e l Sx (k x ,x) = 0, (6) 
where the effective radial potential is obtained by integration over the angles il x 

and m is an arbitrary normalization mass. 



B. Transition matrix element 

The calculation of the transition matrix requires now the overlap of the initial three-body wave function in Eq. (|^) 
and the final state wave function e lfc --R^» u , CT i< J 2 (k x ,x), where k y is the momentum associated with the y coordinate. 
Since k y ■ y = k ■ R we then find 

M^ 2ay (k x ,k y ) = J2 (si<J lS2 a 2 \ Sx a x )M s J ^ Sy<7y (k x ,k y ) (8) 

$x O x 

For simplicity we do not specify the dependence of a (k x , k y ) on the spins si, s 2 and s y of the three particles 

involved in the reaction. 

The analytic form of M^ xS a (fe x , k y ) is 

Mi M aM (k x ,k y ) « 1 J2 E 7 H1 ^ OL K )Y txmix (n k JY iymiy (Q ky ) 

x Y (-iy ,+2S - 2M+e y+ s y- s ^j 2 x j 2 v JLS ( £ jx Jy ) 
^ K ' JXJy \ M —m x —m y J 

x ( h s y \ ( J* ^ x Sx \ 1 ]j \ ^ y 

\ -m y me cr y J \ -m x mt x cr x J | c x y 

I D S x Sy 
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where a = \J2a + 1. In Eq. (||) we have introduced the hyperspherical coordinates in momentum space k 
a K = arctan (k x /k y ), and the numerical function (rt, a re ) ^ s gi ven by 

Kn J 

da(smay* + \cosayy +2 Py hev+i (cos(2a))j t (k yy yf s (k x ,x) 



lkl + kl, 



where v = (K — l x — i y )/2, 2 ' ly + 2 is a Jacobi polynomial, and 



1\ K — 



+ £ x + e y + l)\2(K + 2) 
i> + 4 + f)i> + 4, + §) 



1/2 



(10) 



(11) 



C. Cross sections 



The cross section or momentum distribution is now obtained by squaring the transition matrix, and subsequently 
averaging over initial states and summing up over final states. If we consider the possibility of having a polarized 
beam and measuring the spin projection of the two particles in the final state we can then write 



d 6 a 



M 



mit . (M) ]T W Rn . (a,) WW. (0-2) £ \K^a v (*x, k y 

<T\ G-2 (Ty 



(12) 



where Wmit^M), Wfi n .(ci), and Wfi n .(<72) give the occupation probability of each magnetic substate in the initial and 
final states. When the beam is unpolarizcd and there is no detection of final spin projections these probabilities take 
the values Wa n .(M) = 1/(2J+ 1), Wfin.(ai) = l/(2si + 1), and Wfi n .(o2) = l/(2s2 + 1), and we obtain the unpolarized 
expression for the momentum distributions J19[. 



The volume element dk x dk y is written 



dk x dk y = k x dk x dk x dipk x kydk y dflk y 



(13) 



where k x — k x sin 9k x , k\ — k x cos 9k x , and (9k x , y/c x ) are the polar and azimuthal angles associated to ■ 

The integration of Eq.(|l2|) over £lk y and ipk x can be done analytically. By further numerical integration over k y 

and k x , we get the one-dimensional relative momentum (fcl) distribution of the remaining particles. By integrating 

over k y and k x , we get instead the two-dimensional relative momentum (k x ) distribution. It should be noted that we 
have not specified any coordinate system and the axes x, y, and z are therefore completely arbitrary. Thus, in the 
sudden approximation the longitudinal and transverse momentum distributions are identical. 

After substituting Eqs.(||) and (||) into (12), and integrating analytically over £lk y and ipk x we obtain the expressions 
for the th ree- differential momentum distributions that are collected in appendix A. Three different cases are considered: 

(i) Eq.(Al) describes the process with 100% polarized beam where the polarization of one of the particles in the 
final state is measured, i.e. W init .{M') = S M M', Wn n .(cr'i) = <Wj, Wfi n .(a' 2 ) = l/(2s 2 + 1), 

(ii) Eq.(A2) describes the process with 100% polarized beam without detection of the final spin states, i.e. 



W init .(M') = 6mm>, WW.M) = 1/(2*1 + 1), W Rn .{<j' 2 ) = 1/(2*2 + 1), 

(iii) Eq.(A3) describes the process with unpolarized beam without detection of the final spin states, i.e. Wi n a (M') = 
1/(2.7 + 1), WfinK) = l/(2si + l), W &n .(a' 2 ) = l/(2s 2 + l). 

To compare with the experimental data the momentum distributions should be referred to the center of mass frame 
of the projectile. To do this we construct the momentum p of one of the particles in the final state relative to the 
center of mass of the projectile as a linear combination of k x and k y , i.e. 



p = a,ik x + bik 



(14) 



where i = 1, 2 refer to the two remaining particles. If the Jacobi coordinate x goes from particle 1 to particle 2, then 
a and b take the values 
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fl 171x1712 V /2 x mi f 1 {mi + m 2 )m 3 \ 1/2 
oi = , &i = , 15 

\ TO mi + TO2 / TOl + TO2 \ TO TO-l + TO2 + TO3 / 

/l TO1TO2 \ 1/2 m 2 / 1 (mi +TO 2 )m 3 \ 1/2 
»2 = : , &2 = ; ; : , (16) 

V TO TOl + TO 2 / TOl + TO2 \ TO Vd\ + TO 2 + TO 3 / 

when we compute the momentum distributions of particle 1 and 2, respectively. 
Since the volume element dpdk y — a^dk x dk y we obtain the relation 



d 6 a 1 



dpdky af dk x dk 



(17) 



The only additional complication in this coordinate system is that the numerical function I t y s . («, a„) now depends 



on the relative angle 0^ between p and fe a . Only two angular integrations can then be performed analytically, and 
three integrations must be done numerically in order to get the one-dimensional and two-dimensional momentum 
distributions of one of the particles in the final state relative to the center of mass of the projectile. 

In appendix B we show the analytic expressions for the four-differential momentum distributions relative to the 



center of mass of the projectile for the same three cases as specified in appendix A. They are given by Eqs.(Bl), flB2j) 



and (B3), respectively. 



D. Polarization Observables 



We shall denote the one and two-dimensional momentum distributions when the beam is unpolarized and the spin 
projection of the particles is not measured in the final state by Qi(p") = da/dp" and Qzip^) — da/dp^, respectively. 
The distributions are referre d to the center of mass of the three-body projectile, and obtained by integration over the 
unobserved quantities of Eq.(jB3"|) in appendix B. 

In the same way we shall denote the one and two-dimensional momentum distributions when the beam is 100% 
polarized (with spin projection M) and the spin projection of the particles in the final state is not measured by 
Qf 1 (pjjj = da M /dp" and Q^(p^) = da M /dp ± , respectively. These distributions are obtained by integration of 



Eqs.(B2) in appendix B. 

Finally, <2f f<T (p") = da Mcr /dp" and Q^ 117 {p^) = da Ma /dp 1 - will denote the one and two-dimensional momentum 
distributions for a 100% polarized beam (with spin projection M), and when the spin projectio n of one of the particles 



in the final state is measured to be a. These distributions are obtained by integration of Eqs.(Bl) in appendix B. 
We define now the following asymmetries 



O m - o M ' 

iMM — «f« ~ 



a mm = *x ^ » = 1,2 (18) 



A 



MM' ; 



_ Qf CT - Qf 



M'c 



7=1,2 



(19) 



It is easy to check that the I = I' = term in Eq.(Bl) and the 1 = term in Eq.(B2) are independent of M and a, 
and equal to the totally unpolarized momentum distribution Eq. (|B3| ). This means that the deviations from zero of 
the asymmetries Eqs.(p^|) and (19) are a direct consequence of the polarization of the beam and the measure of the 
spin projection of one of the particles after the fragmentation. 

Direct examination of Eq.(B2) permits us to observe that if only s-waves were involved in the initial and final 
wave functions (£ x = t' x = L x = and £ y = £' y = L v = 0) only the 1=0 term would contribute to the momentum 
distribution. In other words, the asymmetry Eq.(ffq) would vanish. Therefore, the deviation from zero of this 
asymmetry is a signal of non-vanishing components different from s-waves i n t he wave function. This conclusion is 
not true for the asymmetry Eq.(|l9|), since the momentum distribution Eq.(Bl) has additional non- vanishing terms 
even when only s-waves are involved. 
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FIG. 1. The radial wave function fi(p) for n Li when 10 Li has a virtual s-state at 100 keV the lowest-lying p-resonance at 
0.5 MeV. Calculations with 4% of p-wave (solid line), 18% of p-wave (short-dashed line), 26% of p-wave (dotted line), and 32% 
(long-dashed line) in the 10 Li wave function are shown. The inset shows the same wave function when 10 Li has a low-lying 
virtual s-state at 100 keV (dashed line) and when 10 Li has a low-lying p-resonance at 100 keV (solid line). 



III. RESULTS 



In this section we apply the procedure described above to two-neutron halo nuclei, and in particular to n Li. We 
consider neutron removal reactions, and the neutron- 9 Li interaction is include in the final state. We only included 
s- and p-waves in the calculations, since components with higher angular momenta are negligibly small. We use the 
two-body potentials described in ref. pi] , i.e. 

(v c + v s + VtS\2 + V so l nn ■ Snn) exp [— (r/6 n „) 2 ] (20) 



V<$ =V S (1 + 7s s c ■ s„) exp [-(r/b nc ) 2 ] , (21) 
V® = (V t + V£l n c ■ 8«c) exp [-(r/b nc ) 2 ] , (22) 

The nucleoli- nucleon potential V nn is made of gaussians containing central (V c ), spin splitting (V^ s ), spin-orbit 
(V so ) and tensor (Vr) parts where the parameters are fitted to low energy s- and p-wave nucleon-nucleon scattering 
data and s„2 are the spins of the two neutrons, l nn is the two-neutron relative orbital angular momentum, 
Snn = s ni + &n2i and S12 is the usual tensor operator). In the neutron-neutron subsystem more than 98% of the 
probablity is found in s- waves. 

The neutron-core potential V nc is also made of gaussians adjusted to reproduce the experimental binding energy 
and root mean square radius of the three-body projectile (295 ± 35 keV |24| and 3.1 ± 0.3 fm j25|], respectively, for 
11 Li). The different waves are independently fitted, allowing us to adjust separately the energy positions of the virtual 
s-states and the resonances in the neutron-core subsystem. For simplicity the spin splitting term (introduced by the 
parameter 7 S ) has been included only for s- waves in the neutron-core interaction. 
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A. Initial and final wave functions 



The expression we use for the wave function of the three-body halo projectile is given by Eq.(|2|), where the radial 
functions f n (p) and the coefficients C„Ki x e v Ls x s(p) are computed numerically. To illustrate some of the general 
aspects of the three-body wave functions, we show in Fig. |]the /i(p) function for 11 Li. This function dominates the 
expansion in Eq.(|^) as it contains more than 95% of the norm of the total wave function. The neutron-core interaction 
has been fitted to produce a low-lying virtual s-state at 100 keV and the lowest p-resonance at 0.5 MeV. By use of 
the spin-orbit parameter (Vsb~ ) we can move the higher-lying p-resonances up or down and thereby modify the 
p-wave content in the neutron- 9 Li subsystem. Four different situations of varying p-wave content are shown, i.e. 4% 
(solid line), 18% (short-dashed line), 26% (dotted line) and 32% (long-dashed line). We observe that the more p- wave 
content in the neutron- 9 Li subsystem the narrower the radial wave function. This result can be attributed to the 
centrifugal barrier, that keeps the neutrons closer to the core for p-waves. 

When the final state interaction is neglected in the calculation, the momentum distributions are given by the 
square of the Fourier transform of the three-body wave function in Eq. (^) . Therefore a larger p-wave content in the 
neutron-core subsystem will produce a broader momentum distribution. Details about calculations of momentum 
distributions following this method have been presented in Nevertheless, the differences between the three-body 
wave functions in Fig. [I] are not very large and the corresponding momentum distributions would also be very similar. 
In the inset of Fig. |lj we display the same radial three-body wave function j\ (p) for a neutron-core potential producing 
a p- resonance at 100 keV (solid line) and a virtual s-state at 100 keV (dashed line). Again we find that the low- lying 
p- resonance produces a narrower radial wave function, and therefore a broader momentum distribution (without FSI). 
In this case the difference between the wave functions is larger and as a consequence the momentum distributions for 
a low-lying p-resonance and a low-lying virtual s-state will also be clearly distinguishable. 



virtual s-state at 0.1 MeV 




1 2 3 4 5 

*** 

FIG. 2. Continuum radial s-state wave functions in Eq.(Q) when the neutron- 9 Li system has a virtual s-state at 0.1 MeV. 
Three different two-body scattering energies (E2 = h 2 k^/2m) are considered, 0.1 MeV, 1 MeV and 5 MeV. The dashed line 
is the same wave function when the neutron- 9 Li interaction is neglected. The inset show the same wave function when the 
neutron- 9 Li subsystem has a virtual s-state at 1 MeV. 

To analyze the effect of the inclusion of the final state interaction we investigate the behaviour of the distorted 
wave function in Eq.(^). First we consider a neutron-core potential with a virtual s-state at 100 keV, and investigate 
how the corresponding continuum radial s-state wave functions behave for different energies E 2 {= k 2 x h 2 /2m) of the 
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neutron-core system. When the final state interaction is neglected this radial s-state wave function is simply the 
Bessel function jo{k x x), which in Fig. || (the dashed line) is compared to the distorted wave functions (the solid lines) 
for different values of i?2=0.1 MeV, 1 MeV and 5 MeV. At large distances (not shown in the figure) the oscillations 
continue and have the same amplitude. 

The difference from the Bessel function is only significant at small distance and this difference increases in general 
when E2 approaches the energy of the virtual state. (Note, however, the singularity at the origin for vanishing E2-) 
Therefore for small values of k x the overlap in Eq.(Q) is larger than in Eq.Q and since the normalization is the same 
in all cases, the momentum distribution becomes narrower by inclusion of final state interactions. 

The inset in Fig. ^| shows the same plot as in the main part, but for a neutron-core potential with a virtual s-state 
at 1 MeV. The general conclusions are the same, the smaller E2 the larger the effect of the final state interaction. 
However, in this case the difference with the Bessel function is clearly smaller (note the different vertical scales), and 
therefore also the influence of the final state interaction is smaller. We can then conclude that the inclusion of final 
state interaction in s-waves makes the momentum distribution narrower, and the smaller the energy of the virtual 
s-state the larger the effect of the final state interaction. 




1 2 3 4 5 

*** 



FIG. 3. The same as Fig. |2| for p-resonances instead of virtual s-states. 

In Fig. H we show the same quantities as in Fig. |[ but now for p- waves. In the main part of the figure we consider 
a neutron-core potential with a p- resonance at 0.1 MeV. Again the dashed line is the radial part of the p-wave without 
final state interaction (equal to the Bessel function ji(k x x)). When we solve the Schrodinger equation for a p-wave 
and when £/2(= k^h 2 /2m) is far from the resonance energy (we have taken E2 = 5 MeV) the difference with the 
Bessel function is very small, and the effect of the final state interaction is almost negligible. When E 2 approaches 
the energy of the resonance the p-wave becomes more and more localized, in such a way that it gets a maximum for 
E2 equal to the energy of the resonance (0.1 MeV in this case). As a consequence, the overlap in Eq.(^) is maximum 

for k x — \J 2mE TCS /h 2 , and the momentum distributions arc enhanced around that value of k x . 

The inset of Fig. ^ shows the same as in the main part for E Tes = 1 MeV. The behaviour of the p-wave is identical to 
the previous one. The wave is very much localized and gets its maximum value for E 2 = E rcs . However, the maximum 
is now clearly smaller than in the previous case of E les = 0.1 MeV, and therefore also the effect on the momentum 
distributions is less pronounced. Therefore, the inclusion of the final state interaction in the p-wave will produce 

an increase in the momentum distributions around k x = \J 2mE Tes /h 2 . This increase is especially important when 
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the neutron-core potential has a low-lying p-resonance, and again, due to the normalization, makes the momentum 
distributions narrower. 

Summarizing, the inclusion of final state interaction highly modifies the momentum distributions for s- and p- waves. 
This change is clearly more significant for low-lying virtual s-states and low-lying p-resonances. For s-waves the 

distributions arc enhanced at the origin (k x = 0), and for p- waves the enhancement is produced at k x = \J 2mE res /h 2 . 
In both cases the momentum distributions become narrower by inclusion of the final state interaction. 

Before concluding this subsection we must emphasize that the momentum distributions discussed so far were 
expressed in terms of k x , i.e. the relative momentum between the two remaining particles. From Eqs.(|l4]) - ( |l6| ) 
we see that the momentum of the neutron relative to the center of mass of the projectile is very close to k x (the b 
coefficient is small due to the large mass of the core). Therefore the above conclusions about the effect of the final 
state interaction roughly remain unchanged for the neutron momentum distributions relative to the center of mass of 
the projectile. 

For the core momentum distribution the hi coefficient in Eq.([l4|) is large, and the effect of the final state interaction 
is much smaller when we refer momentum distributions of the core to the center of mass of the projectile. The reason 
is simply that the center of mass of the core is fairly close to the center of mass of the total three-body system. 

B. Asymmetries 

Experimental data suggest that the neutron unbound 10 Li system has a p-resonance around 0.5 MeV |^6| and 
a low- lying (< 0.2 MeV) virtual s-state or p-resonance p] , p6| |28| . In all the computations reported in this work 
the neutron-core potential is fitted to produce a p-resonance at 0.5 MeV in 10 Li. For the lowest-lying state we 
investigate several possibilities in order to analyze the dependence of the polarization observables on the properties 
of the lowest-lying state. 
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FIG. 4. Asymmetry A t (p") for neutron momentum distributions in a neutron removal Li fragmentation reaction. A 
low-lying state is assumed in 10 Li at 50 keV (solid line), 100 keV (dashed line) or at 150 keV (dot-dashed line). The cases of 
low-lying p-resonances and low-lying virtual s-states are both shown (see figure) . 

In Fig. ^ we show the asymmetry A^ 2 1 ^ 2 (p") for neutron momentum distributions after a neutron removal n Li 
fragmentation reaction. A low-lying state (s- or p-wave) is present in 10 Li at 50 keV (solid line), 100 keV (dashed 
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line) or at 150 keV (dot-dashed line). For low-lying p-resonances we observe for all the three energies two peaks in the 
asymmetry. The first peak at low momentum corresponds to the lowest-lying p-resonance. The second peak is in all 
the three cases generated by the p-resonance at 0.5 MeV. The connection between the position of the differe nt min ima 
and maxima of the asymmetry and the energy of the p-resonance may be seen from the relations k x = p r ^Jm/ \x and 

k x = \J 2mE Tes /h 2 , where p r is the relative neutron-core momentum, E Ies is the energy of the resonance and fi is the 
reduced mass of the neutron-core subsystem. The different p-resonance energies in the figure correspond to values p r 
of 10 MeV/c, 14 MeV/c, and 17 MeV/c, respectively and an energy of 0.5 MeV corresponds to p r « 30 MeV/c. 

As discussed in section 2, the momentum distributions are enhanced around the resonance value of k x , that is very 
close to the momentum of the neutron relative to the center of mass of the three-body projectile. This increase in 
the momentum distributions produces the oscillations in the asymmetry. As a consequence, when the energy of the 
resonance increases the peak in the asymmetry is displaced towards higher momenta. Since we plot the asymmetries 
as functions of the longitudinal component of the neutron momentum relative to the center of mass of the projectile 
(that is not exactly p r , but very close), the maxima and minima of the asymmetry do not exactly coincide with these 
numbers. However, the connection to the energy of the resonances in the neutron-core subsystem is very clear. 

Inspecting the asymmetry in Fig. |^ where a low-lying virtual s-state is present in 10 Li we observe that the first peak 
completely disappears. This fact was expected, since in section 2 we showed that the asymmetry Eq.([l8|) vanishes 

when only s- waves are involved in the wave functions. In other words, the asymmetry A^ 2 ^ 2 (p") does not feel 
the presence of virtual s-states in the neutron-core subsystem, and it is only sensitive to the structure of the non 
zero angular momnetum resonances. Since only s-waves are changed for the three different cases considered, the 
asymmetries are essentially indistinguishible. 




FIG. 5. The same as in Fig. ^ for the asymmetry A^ 2 1/ ' 2 ' 3/ ' 2 3,72 (p"). 

Let us now consider the case where the spin projection of one of the particles in the final state is measured. 
In particular we assume that the spin projection of the core is measured to be 3/2, and construct the asymmetry 

3/2 l/2 - 3/2 3/2 II n 

A ± (p") for neutron momentum distributions in a neutron removal reaction. The result is shown in Fig. a 

for the cases corresponding to Fig. ||. Again, for low-lying p-resonances there are two peaks in the asymmetry, one of 
them connected to the energy of the low- lying resonance and the other peak is related to the resonance at 0.5 MeV. 

For low-lying virtual s-states the first peak disappears, and only the peak associated with the p-resonance at 0.5 
MeV remains. Now the s-waves also contribute to the asymmetry, but the contribution is localized in the region close 
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to the origin where the curves now are different for the three energies of the low-lying virtual s-state. Therefore the 
asymmetries shown in Fig. [| display the same features as those of Fig. ||. The main difference is that the amplitude 
of the oscillations of the asymmetries clearly is larger in the asymmetries shown in Fig. 0. 




3/21/211 11 

FIG. 6. Asymmetry A t (p") for neutron momentum distributions in a neutron removal Li fragmentation reaction. A 
low-lying virtual s-state and a p-resonance are assumed in 10 Li at 100 keV and 0.5 MeV, respectively. The p-wave content in 
10 Li is considered to be 4% (solid line), 18% (dashed line), and 32% (dot-dashed line). 



We have now concluded that the asymmetries are sensitive essentially only to the non-zero angular momentum 
partial waves in the initial and final state wave function. Thus we also expect sensitivity to the content of the 
different waves in the two-body subsystem. In particular, for n Li the asymmetries should depend on the p-wave 
content in 10 Li subsystem (the neutron-neutron subsystem contains in any case almost exclusively s-waves). In Fig. 

U we show the same asymmetry, A^ 2 ^ 2 (p"), as in Fig. || but now assuming a virtual s-state at 100 keV and a 
p-resonance at 0.5 MeV in 10 Li. By use of the spin-orbit parameter Vj^ 1 in the neutron-core potential in Eq.(^2|) we 
modify the p-wave content in the neutron- 9 Li subsystem. 

The solid line, dashed line, and dot-dashed line are the resulting calculations with 4%, 18%, and 32% of p-wave 
content in the 10 Li wave function. In the three cases we observe the peak in the asymmetry around 30 MeV/c 
corresponding to the p-resonance at 0.5 MeV. The second peak in the asymmetry corresponds to the next p-resonance, 
closer to the first peak for a larger p-wave content, since the increase in the p-wave content is obtained by decreasing 
the energies of the p-resonances. 

We observe that the lower the p-wave content the smaller the amplitudes of the oscillations of the asymmetry, see 
the solid line corresponding to the p-wave content of 4%. This is related to the fact that the asymmetry must vanish 
when the p-waves are totally absent, However, it should be mentioned that the denominator of the asymmetry in 
Eq.(|l8|) is the unpolarized momentum distribution, that is known to be narrower for lower p-wave contents ]l9| . It is 
then possible to find oscillations with similar amplitude for different values of the p-wave content. 
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FIG. 7. The same as in Fig. ^ for the asymmetry A^ 2 1//2,3//2 3//2 (p"). 

3/2 l/2'3/2 3/2 II C^l 
we show the same asymmetry, A x (p 11 ), as in Fig. |6|, i.e. one-dimensional neutron momentum 

distributions in a neutron removal reaction. As in Fig. |^ the spin projection of the core is assumed to be measured 

to have a value of 3/2 in the final state. Again, in this case the asymmetry shows essentially the same features as in 

Fig. |^, but with the advantage of larger amplitudes. The values of the asymmetries can vary now between 0.8 and 

more than 1.3. 
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FIG. 9. The same as in Fig. | for the asymmetry A\ /2 1/2;3/2 3/2 (p^). 

All the asymmetries we have discussed so far are obtained from one-dimensional neutron momentum distributions. 
The same conclusions can be obtained from the corresponding two-dimensional neutron momentum distributions 
shown for completeness in Figs. ||, 0, |l^, O. First we note that when only the beam is polarized the asymmetries 
are not sensitive to the structure of the virtual s-states in the two-body subsystems. When only the energy of the 
s-states is modified in the 10 Li wave function the asymmetries do not change. Second, the minima and maxima of 
the asymmetries provide information about the p-resonance structure in the 10 Li two-body subsystem. Third, the 
asymmetries are clearly dependent on the p-wave content of the 10 Li wave function. Fourth, the asymmetries obtained 
assuming that the spin projection of the core has been measured in the final state (3/2 in our case) show a clearly 
larger effect than those obtained assuming that only the beam is polarized. 

In principle the same asymmetries could also be obtained from core momentum distributions and from neutron 
momentum distributions after core break-up. However these two cases are not especially interesting. The former 
because core momentum distributions are very little affected by the final state interaction, and then the asymmetries 
do not show very much structure. The latter because after core break-up the two neutrons survive in the final state, 
and the three-body wave function written in the Jacobi set where the x coordinate is drawn between the two neutrons 
contains almost exclusively s-waves, to which the asymmetries are not sensitive. 

IV. SUMMARY AND CONCLUSIONS 

We have formulated a detailed model to investigate momentum distributions of particles emerging after nuclear 
break-up reactions of fast polarized three-body halo nuclei. The three main ingredients are: 

(i) The three-body description of the initial nuclear wave function. This assumes that the particles or clusters are 
inert throughout the process and only distances larger than the cluster radii can be well described. Such models have 
been successfully applied to various light nuclei. 

(ii) The sudden approximation where the target nucleus instantaneously removes one of the particles from the three- 
body projectile without disturbing the other two. This approximation is very well suited for short range interactions 
and fast processes where the beam velocity is much larger than the velocity of the intrinsic motion. The radioactive 
beams are available with high energies exceedingly well fulfilling this criterium. 

(iii) The final state relative wave function for the two remaining particles must include effects of their mutual 
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interaction. The resulting distorted wave function may differ substantially from the plane wave solution when no two- 
body interaction is present. These effects are crucial when virtual s-states or higher angular momentum resonances 
occur at energies below about 1 MeV in the final state two-body system. 
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It is rather difficult to construct a theoretical system where assumption (i) is valid without the necessity of using 
(iii). The reason is that the three-body structure only is appropriate when the the particles are fairly loosely bound, 
since they otherwise would overlap and destroy the cluster property. As a consequence the two-body subsystems must 
have fairly weak effective attractions yet they must be sufficiently strong to create a three-body bound state when 
working together. The result is low-lying continuum two-body structure. 

The same two-body interactions are therefore responsible for both the initial three-body bound state and for the 
distorted final state wave function. A consistent treatment of initial and final stages of the process is therefore essential 
and our model meets this crucial requirement. In other words the two-body continuum structure is decisive for both 
the three-body halo structure and the momentum distributions after fragmentation. 

We first describe the method to compute the momentum distributions in nuclear break-up processes. The overlap 
between initial and final state wave function is the essential quantity. The initial bound state three-body halo wave 
function is obtained by solving the Faddeev equations by use of an adiabatic hyperspherical expansion technique. 
This method has been tested previously and shown to be very efficient for ground states or low-lying excited states 
with small two-body relative orbital angular momenta. The final state distorted two-body continuum wave function 
is expanded in partial waves and the related radial wave functions are calculated numerically. They are normalized 
at large distances to coincide with the plane wave solutions. The resulting overlap matrix with different angular 
momentum projections is then available. 

The cross sections or the observable momentum distributions are then essentially only the absolute square of the 
overlap integrated over unobserved quantities. These expressions contain numerically computed quantities. They 
appear to be rather complicated especially due to the use of the coordinate system where the center of mass of the 
three-body projectile is at rest and the dedicated goal of computing momentum distributions for polarized projectiles. 

First we specify the two-body interactions appropriate for the 11 Li system considered as two neutrons surrounding 
the 9 Li-core. The neutron-neutron interaction is parametrized to reproduce low-energy scattering properties and then 
maintained throughout the present investigation. The neutron-core interaction is adjusted to reproduce the measured 
11 Li binding energy and root mean square radius. The remaining free parameters are used to place virtual s-states and 
p-resonances at various energies close to the threshold resulting in different 11 Li-structures and different momentum 
distributions. 

The final state two-body s-wave functions all have maxima at the origin. The highest value occurs when the 
scattering energy equals to the energy of the virtual s-state. The state is therefore peaked around zero and this 
localization is strongest around the energy of the virtual state. The p-state two-body wave function is peaked at an 
energy corresponding to the p-resonance energy. Without resonances and virtual states the wave functions are simply 
Bessel functions with the related oscillations and without these pronounced peak structures. 

Neutron momentum distributions after neutron removal are computed for polarized 11 Li and when the polarization 
of the final two-body state is measured. Several polarization observables are constructed and expressed in terms of 
quantities, called asymmetries, deviating from zero only when the momentum distributions depend on the polarization. 
Several continuum structures of 10 Li are considered, but in all cases we maintain the binding energy and root mean 
square radius of 11 Li. 

The inclusion of the final state interaction is an essential ingredient in the computation, determining the behaviour 
of the asymmetries. In fact, when they are not included the structure in the polarization observables completely 
disappears. The core momentum distributions are little affected by the final state interaction, and therefore the 
asymmetries obtained from them exhibit very little structure. 

From the asymmetries constructed from neutron momentum distributions we have observed the following features: 
If only the beam is polarized the asymmetries vanish when only s-waves contribute. This makes the asymmetries 
obtained from a core break-up reaction from 11 Li not relevant, because almost exclusively s-waves are involved in this 
fragmentation process. Much more interesting are the asymmetries obtained from neutron momentum distributions 
after neutron removal reactions. The rest of the waves plays an important role, and the asymmetries become a perfect 
tool to investigate the I > O-resonance structure of the neutron-core subsystem. These observables then contain 
information about the continuum structure of the two-body subsystems. In fact low-lying resonances (p-resonances 
in the 11 Li case) in the neutron-core system show up as an extremum at the momentum value corresponding to the 
resonance energy. An additonal measurement of the two-body final state polarization enhances the amplitudes of the 
asymmetries, but otherwise leaves essentially unchanged the structure of the asymmetry. 

APPENDIX A: RELATIVE MOMENTUM DISTRIBUTIONS 

In this appendix we give the analytic expressions for the relative momentum distributions of the two particles in 
the final state. We consider three different cases: 
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(i) 100% polarized beam and spin projection of one of the particles in the final state measured: 

5^ £<-»'-'" (i -I, {)»■*-*■> 
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(ii) 100% polarized beam and no spins projections measured in the final state: 

This amounts to multiply the previous expression by l/(2s\ + 1) and sum over u\. The main consequence of the 
summation is that the index V is zero, and then L x = I. Then we get the expression 
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(hi) Unpolarized beam and no spins projections measured in the final state: 

Multiplying by 1/(2 J + 1) and summing over M we get the totally unpolarized momentum distribution. From the 
summation we obtain that 1 = 0, and the analytic expression is 
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APPENDIX B: MOMENTUM DISTRIBUTIONS RELATIVE TO THE CENTER OF MASS OF THE 

PROJECTILE 

In this appendix we give the analytic expressions for the momentum distributions of the two particles in the final 
state relative to the center of mass of the incident three-body system. The same three cases as in the previous 
appendix are considered: 

(i) 100% polarized beam and spin projection of one of the particles in the final state measured: 
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(ii) 100% polarized beam and no spins projections measured in the final state: 



dkyd6' k dp^dpW 



|j^(-i) J - M / 2 Px(co S e p ) 



j j i 

M -M 



E E E ^ v ^'^l*j*Cjy 2 PLWyLUSS< 

l*Sx3*LS l^j'^L'S' l y j y l' y3 ' y 



^x s x3 x 



J'L'S' 





jx 


jy) r j 


Jx Jy I 


is 


4 




£' £' \ 




Sx 


Sy) {S' 


Sx Sy J 



L X Ly 



£ x ^X. L-y 



'X ^X ^X 





f P' T 

l v v y 




jy jy 



jx j x L x 
jy jy L y 
J J I 



jx jx L x 



E 



L X Ly I 



(L*-w v+ Nxy. p ^ (coseL) P -^ (cos(? , 



N x -N x J V (L x + N x y.(L y — N x )\ L * 



(B2) 



(iii) Unpolarized beam and no spins projections measured in the final state 
d 4 a 
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All the three previous expressions are valid when a < 0. For a > the summation over N x should include an extra 
factor (-1)^. 

If we assume that one of the particles in the final state, say particle 1, has infinite mass, then the momentum 
distribution of particle 2 relative to the t hree- b ody center of mass should coincide with the momentum distrib utio n 
relative to particle 1. In other words, Eqs.(Al), (A2), and (A3) should be recovered by integrating over 0' k Eqs.(Bl), 

(B2), and (|B3|), respectively. In this case 6, = 0, see Eq. (|l6|) , p = ciik x and therefore 8' k = and P^"(l) = <5jv x ,o- 
Now k and a K are independent of 9' k and the integration over 9' k can easily be done analytically leading directly to 
Eqs.Q, (|A§), and Q. 
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